home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
GameStar 2004 April
/
Gamestar_61_2004-04_dvdb.iso
/
DVDStar
/
Editace
/
hltp.exe
/
{app}
/
Source Code
/
Zoners Half-Life Tools
/
hlrad
/
vismatrixutil.cpp
< prev
Wrap
C/C++ Source or Header
|
2002-12-09
|
22KB
|
713 lines
#include "qrad.h"
funcCheckVisBit g_CheckVisBit = NULL;
unsigned g_total_transfer = 0;
unsigned g_transfer_index_bytes = 0;
unsigned g_transfer_data_bytes = 0;
#define COMPRESSED_TRANSFERS
//#undef COMPRESSED_TRANSFERS
int FindTransferOffsetPatchnum(transfer_index_t* tIndex, const patch_t* const patch, const unsigned patchnum)
{
//
// binary search for match
//
int low = 0;
int high = patch->iIndex - 1;
int offset;
while (1)
{
offset = (low + high) / 2;
if ((tIndex[offset].index + tIndex[offset].size) < patchnum)
{
low = offset + 1;
}
else if (tIndex[offset].index > patchnum)
{
high = offset - 1;
}
else
{
unsigned x;
unsigned int rval = 0;
transfer_index_t* pIndex = tIndex;
for (x = 0; x < offset; x++, pIndex++)
{
rval += pIndex->size + 1;
}
rval += patchnum - tIndex[offset].index;
return rval;
}
if (low > high)
{
return -1;
}
}
}
#ifdef COMPRESSED_TRANSFERS
static unsigned GetLengthOfRun(const transfer_raw_index_t* raw, const transfer_raw_index_t* const end)
{
unsigned run_size = 0;
while (raw < end)
{
if (((*raw) + 1) == (*(raw + 1)))
{
raw++;
run_size++;
if (run_size >= MAX_COMPRESSED_TRANSFER_INDEX_SIZE)
{
return run_size;
}
}
else
{
return run_size;
}
}
return run_size;
}
static transfer_index_t* CompressTransferIndicies(transfer_raw_index_t* tRaw, const unsigned rawSize, unsigned* iSize)
{
unsigned x;
unsigned size = rawSize;
unsigned compressed_count = 0;
transfer_raw_index_t* raw = tRaw;
transfer_raw_index_t* end = tRaw + rawSize - 1; // -1 since we are comparing current with next and get errors when bumping into the 'end'
transfer_index_t CompressedArray[MAX_PATCHES]; // somewhat big stack object (1 Mb with 256k patches)
transfer_index_t* compressed = CompressedArray;
for (x = 0; x < size; x++, raw++, compressed++)
{
compressed->index = (*raw);
compressed->size = GetLengthOfRun(raw, end); // Zero based (count 0 still implies 1 item in the list, so 256 max entries result)
raw += compressed->size;
x += compressed->size;
compressed_count++; // number of entries in compressed table
}
*iSize = compressed_count;
if (compressed_count)
{
unsigned compressed_array_size = sizeof(transfer_index_t) * compressed_count;
transfer_index_t* rval = (transfer_index_t*)AllocBlock(compressed_array_size);
ThreadLock();
g_transfer_index_bytes += compressed_array_size;
ThreadUnlock();
memcpy(rval, CompressedArray, compressed_array_size);
return rval;
}
else
{
return NULL;
}
}
#else
static transfer_index_t* CompressTransferIndicies(const transfer_raw_index_t* tRaw, const unsigned rawSize, unsigned* iSize)
{
unsigned x;
unsigned size = rawSize;
unsigned compressed_count = 0;
transfer_raw_index_t* raw = tRaw;
transfer_raw_index_t* end = tRaw + rawSize;
transfer_index_t CompressedArray[MAX_PATCHES]; // somewhat big stack object (1 Mb with 256k patches)
transfer_index_t* compressed = CompressedArray;
for (x = 0; x < size; x++, raw++, compressed++)
{
compressed->index = (*raw);
compressed->size = 0;
compressed_count++; // number of entries in compressed table
}
*iSize = compressed_count;
if (compressed_count)
{
unsigned compressed_array_size = sizeof(transfer_index_t) * compressed_count;
transfer_index_t* rval = AllocBlock(compressed_array_size);
ThreadLock();
g_transfer_index_bytes += compressed_array_size;
ThreadUnlock();
memcpy(rval, CompressedArray, compressed_array_size);
return rval;
}
else
{
return NULL;
}
}
#endif
/*
* =============
* MakeScales
*
* This is the primary time sink.
* It can be run multi threaded.
* =============
*/
#ifdef SYSTEM_WIN32
#pragma warning(push)
#pragma warning(disable: 4100) // unreferenced formal parameter
#endif
void MakeScales(const int threadnum)
{
int i;
unsigned j;
vec3_t delta;
vec_t dist;
int count;
float trans;
patch_t* patch;
patch_t* patch2;
float send;
vec3_t origin;
vec_t area;
const vec_t* normal1;
const vec_t* normal2;
vec_t total;
transfer_raw_index_t* tIndex;
transfer_data_t* tData;
transfer_raw_index_t* tIndex_All = (transfer_raw_index_t*)AllocBlock(sizeof(transfer_index_t) * MAX_PATCHES);
transfer_data_t* tData_All = (transfer_data_t*)AllocBlock(sizeof(transfer_data_t) * MAX_PATCHES);
count = 0;
while (1)
{
i = GetThreadWork();
if (i == -1)
break;
patch = g_patches + i;
patch->iIndex = 0;
patch->iData = 0;
total = 0.0;
tIndex = tIndex_All;
tData = tData_All;
VectorCopy(patch->origin, origin);
normal1 = getPlaneFromFaceNumber(patch->faceNumber)->normal;
area = patch->area;
// find out which patch2's will collect light
// from patch
for (j = 0, patch2 = g_patches; j < g_num_patches; j++, patch2++)
{
vec_t dot1;
vec_t dot2;
#ifdef HLRAD_HULLU
vec3_t transparency = {1.0,1.0,1.0};
#endif
#ifdef HLRAD_HULLU
if (!g_CheckVisBit(i, j, transparency) || (i == j))
#else
if (!g_CheckVisBit(i, j) || (i == j))
#endif
{
continue;
}
normal2 = getPlaneFromFaceNumber(patch2->faceNumber)->normal;
// calculate transferemnce
VectorSubtract(patch2->origin, origin, delta);
dist = VectorNormalize(delta);
dot1 = DotProduct(delta, normal1);
dot2 = -DotProduct(delta, normal2);
trans = (dot1 * dot2) / (dist * dist); // Inverse square falloff factoring angle between patch normals
#ifdef HLRAD_HULLU
trans = trans * VectorAvg(transparency); //hullu: add transparency effect
#endif
if (trans >= 0)
{
send = trans * patch2->area;
// Caps light from getting weird
if (send > 0.4f)
{
trans = 0.4f / patch2->area;
send = 0.4f;
}
total += send;
// scale to 16 bit (black magic)
trans = trans * area * INVERSE_TRANSFER_SCALE;
if (trans >= TRANSFER_SCALE_MAX)
{
trans = TRANSFER_SCALE_MAX;
}
}
else
{
#if 0
Warning("transfer < 0 (%f): dist=(%f)\n"
" dot1=(%f) patch@(%4.3f %4.3f %4.3f) normal(%4.3f %4.3f %4.3f)\n"
" dot2=(%f) patch@(%4.3f %4.3f %4.3f) normal(%4.3f %4.3f %4.3f)\n",
trans, dist,
dot1, patch->origin[0], patch->origin[1], patch->origin[2], patch->normal[0], patch->normal[1],
patch->normal[2], dot2, patch2->origin[0], patch2->origin[1], patch2->origin[2],
patch2->normal[0], patch2->normal[1], patch2->normal[2]);
#endif
trans = 0.0;
}
*tData = trans;
*tIndex = j;
tData++;
tIndex++;
patch->iData++;
count++;
}
// copy the transfers out
if (patch->iData)
{
unsigned data_size = patch->iData * sizeof(transfer_data_t);
patch->tData = (transfer_data_t*)AllocBlock(data_size);
patch->tIndex = CompressTransferIndicies(tIndex_All, patch->iData, &patch->iIndex);
hlassume(patch->tData != NULL, assume_NoMemory);
hlassume(patch->tIndex != NULL, assume_NoMemory);
ThreadLock();
g_transfer_data_bytes += data_size;
ThreadUnlock();
//
// normalize all transfers so exactly 50% of the light
// is transfered to the surroundings
//
total = 0.5 / total;
{
unsigned x;
transfer_data_t* t1 = patch->tData;
transfer_data_t* t2 = tData_All;
for (x = 0; x < patch->iData; x++, t1++, t2++)
{
(*t1) = (*t2) * total;
}
}
}
}
FreeBlock(tIndex_All);
FreeBlock(tData_All);
ThreadLock();
g_total_transfer += count;
ThreadUnlock();
}
#ifdef SYSTEM_WIN32
#pragma warning(pop)
#endif
/*
* =============
* SwapTransfersTask
*
* Change transfers from light sent out to light collected in.
* In an ideal world, they would be exactly symetrical, but
* because the form factors are only aproximated, then normalized,
* they will actually be rather different.
* =============
*/
void SwapTransfers(const int patchnum)
{
patch_t* patch = &g_patches[patchnum];
transfer_index_t* tIndex = patch->tIndex;
transfer_data_t* tData = patch->tData;
unsigned x;
for (x = 0; x < patch->iIndex; x++, tIndex++)
{
unsigned size = (tIndex->size + 1);
unsigned patchnum2 = tIndex->index;
unsigned y;
for (y = 0; y < size; y++, tData++, patchnum2++)
{
patch_t* patch2 = &g_patches[patchnum2];
if (patchnum2 > patchnum)
{ // done with this list
return;
}
else if (!patch2->iData)
{ // Set to zero in this impossible case
Log("patch2 has no iData\n");
(*tData) = 0;
continue;
}
else
{
transfer_index_t* tIndex2 = patch2->tIndex;
transfer_data_t* tData2 = patch2->tData;
int offset = FindTransferOffsetPatchnum(tIndex2, patch2, patchnum);
if (offset >= 0)
{
transfer_data_t tmp = *tData;
*tData = tData2[offset];
tData2[offset] = tmp;
}
else
{ // Set to zero in this impossible case
Log("FindTransferOffsetPatchnum returned -1 looking for patch %d in patch %d's transfer lists\n",
patchnum, patchnum2);
(*tData) = 0;
return;
}
}
}
}
}
#ifdef HLRAD_HULLU
/*
* =============
* MakeScales
*
* This is the primary time sink.
* It can be run multi threaded.
* =============
*/
#ifdef SYSTEM_WIN32
#pragma warning(push)
#pragma warning(disable: 4100) // unreferenced formal parameter
#endif
void MakeRGBScales(const int threadnum)
{
int i;
unsigned j;
vec3_t delta;
vec_t dist;
int count;
float trans[3];
float trans_one;
patch_t* patch;
patch_t* patch2;
float send;
vec3_t origin;
vec_t area;
const vec_t* normal1;
const vec_t* normal2;
vec_t total;
transfer_raw_index_t* tIndex;
rgb_transfer_data_t* tRGBData;
transfer_raw_index_t* tIndex_All = (transfer_raw_index_t*)AllocBlock(sizeof(transfer_index_t) * MAX_PATCHES);
rgb_transfer_data_t* tRGBData_All = (rgb_transfer_data_t*)AllocBlock(sizeof(rgb_transfer_data_t) * MAX_PATCHES);
count = 0;
while (1)
{
i = GetThreadWork();
if (i == -1)
break;
patch = g_patches + i;
patch->iIndex = 0;
patch->iData = 0;
total = 0.0;
tIndex = tIndex_All;
tRGBData = tRGBData_All;
VectorCopy(patch->origin, origin);
normal1 = getPlaneFromFaceNumber(patch->faceNumber)->normal;
area = patch->area;
// find out which patch2's will collect light
// from patch
for (j = 0, patch2 = g_patches; j < g_num_patches; j++, patch2++)
{
vec_t dot1;
vec_t dot2;
vec3_t transparency = {1.0,1.0,1.0};
if (!g_CheckVisBit(i, j, transparency) || (i == j))
{
continue;
}
normal2 = getPlaneFromFaceNumber(patch2->faceNumber)->normal;
// calculate transferemnce
VectorSubtract(patch2->origin, origin, delta);
dist = VectorNormalize(delta);
dot1 = DotProduct(delta, normal1);
dot2 = -DotProduct(delta, normal2);
trans_one = (dot1 * dot2) / (dist * dist); // Inverse square falloff factoring angle between patch normals
VectorFill(trans, trans_one);
VectorMultiply(trans, transparency, trans); //hullu: add transparency effect
if (VectorAvg(trans) >= 0)
{
/////////////////////////////////////////RED
send = trans[0] * patch2->area;
// Caps light from getting weird
if (send > 0.4f)
{
trans[0] = 0.4f / patch2->area;
send = 0.4f;
}
total += send / 3.0f;
/////////////////////////////////////////GREEN
send = trans[1] * patch2->area;
// Caps light from getting weird
if (send > 0.4f)
{
trans[1] = 0.4f / patch2->area;
send = 0.4f;
}
total += send / 3.0f;
/////////////////////////////////////////BLUE
send = trans[2] * patch2->area;
// Caps light from getting weird
if (send > 0.4f)
{
trans[2] = 0.4f / patch2->area;
send = 0.4f;
}
total += send / 3.0f;
// scale to 16 bit (black magic)
VectorScale(trans, area * INVERSE_TRANSFER_SCALE, trans);
if (trans[0] >= TRANSFER_SCALE_MAX)
{
trans[0] = TRANSFER_SCALE_MAX;
}
if (trans[1] >= TRANSFER_SCALE_MAX)
{
trans[1] = TRANSFER_SCALE_MAX;
}
if (trans[2] >= TRANSFER_SCALE_MAX)
{
trans[2] = TRANSFER_SCALE_MAX;
}
}
else
{
#if 0
Warning("transfer < 0 (%4.3f %4.3f %4.3f): dist=(%f)\n"
" dot1=(%f) patch@(%4.3f %4.3f %4.3f) normal(%4.3f %4.3f %4.3f)\n"
" dot2=(%f) patch@(%4.3f %4.3f %4.3f) normal(%4.3f %4.3f %4.3f)\n",
trans[0], trans[1], trans[2], dist,
dot1, patch->origin[0], patch->origin[1], patch->origin[2], patch->normal[0], patch->normal[1],
patch->normal[2], dot2, patch2->origin[0], patch2->origin[1], patch2->origin[2],
patch2->normal[0], patch2->normal[1], patch2->normal[2]);
#endif
VectorFill(trans,0.0);
}
VectorCopy(trans, *tRGBData);
*tIndex = j;
tRGBData++;
tIndex++;
patch->iData++;
count++;
}
// copy the transfers out
if (patch->iData)
{
unsigned data_size = patch->iData * sizeof(rgb_transfer_data_t);
patch->tRGBData = (rgb_transfer_data_t*)AllocBlock(data_size);
patch->tIndex = CompressTransferIndicies(tIndex_All, patch->iData, &patch->iIndex);
hlassume(patch->tRGBData != NULL, assume_NoMemory);
hlassume(patch->tIndex != NULL, assume_NoMemory);
ThreadLock();
g_transfer_data_bytes += data_size;
ThreadUnlock();
//
// normalize all transfers so exactly 50% of the light
// is transfered to the surroundings
//
total = 0.5 / total;
{
unsigned x;
rgb_transfer_data_t* t1 = patch->tRGBData;
rgb_transfer_data_t* t2 = tRGBData_All;
for (x = 0; x < patch->iData; x++, t1++, t2++)
{
VectorScale( *t2, total, *t1 );
}
}
}
}
FreeBlock(tIndex_All);
FreeBlock(tRGBData_All);
ThreadLock();
g_total_transfer += count;
ThreadUnlock();
}
#ifdef SYSTEM_WIN32
#pragma warning(pop)
#endif
/*
* =============
* SwapTransfersTask
*
* Change transfers from light sent out to light collected in.
* In an ideal world, they would be exactly symetrical, but
* because the form factors are only aproximated, then normalized,
* they will actually be rather different.
* =============
*/
void SwapRGBTransfers(const int patchnum)
{
patch_t* patch = &g_patches[patchnum];
transfer_index_t* tIndex = patch->tIndex;
rgb_transfer_data_t* tRGBData= patch->tRGBData;
unsigned x;
for (x = 0; x < patch->iIndex; x++, tIndex++)
{
unsigned size = (tIndex->size + 1);
unsigned patchnum2 = tIndex->index;
unsigned y;
for (y = 0; y < size; y++, tRGBData++, patchnum2++)
{
patch_t* patch2 = &g_patches[patchnum2];
if (patchnum2 > patchnum)
{ // done with this list
return;
}
else if (!patch2->iData)
{ // Set to zero in this impossible case
Log("patch2 has no iData\n");
VectorFill(*tRGBData, 0);
continue;
}
else
{
transfer_index_t* tIndex2 = patch2->tIndex;
rgb_transfer_data_t* tRGBData2 = patch2->tRGBData;
int offset = FindTransferOffsetPatchnum(tIndex2, patch2, patchnum);
if (offset >= 0)
{
rgb_transfer_data_t tmp;
VectorCopy(*tRGBData, tmp)
VectorCopy(tRGBData2[offset], *tRGBData);
VectorCopy(tmp, tRGBData2[offset]);
}
else
{ // Set to zero in this impossible case
Log("FindTransferOffsetPatchnum returned -1 looking for patch %d in patch %d's transfer lists\n",
patchnum, patchnum2);
VectorFill(*tRGBData, 0);
return;
}
}
}
}
}
#endif /*HLRAD_HULLU*/
#ifndef HLRAD_HULLU
void DumpTransfersMemoryUsage()
{
Log("Transfer Lists : %u transfers\n Indices : %u bytes\n Data : %u bytes\n",
g_total_transfer, g_transfer_index_bytes, g_transfer_data_bytes);
}
#else
//More human readable numbers
void DumpTransfersMemoryUsage()
{
if(g_total_transfer > 1000*1000)
Log("Transfer Lists : %11u : %7.2fM transfers\n", g_total_transfer, g_total_transfer/(1000.0f*1000.0f));
else if(g_total_transfer > 1000)
Log("Transfer Lists : %11u : %7.2fk transfers\n", g_total_transfer, g_total_transfer/1000.0f);
else
Log("Transfer Lists : %11u transfers\n", g_total_transfer);
if(g_transfer_index_bytes > 1024*1024)
Log(" Indices : %11u : %7.2fM bytes\n", g_transfer_index_bytes, g_transfer_index_bytes/(1024.0f * 1024.0f));
else if(g_transfer_index_bytes > 1024)
Log(" Indices : %11u : %7.2fk bytes\n", g_transfer_index_bytes, g_transfer_index_bytes/1024.0f);
else
Log(" Indices : %11u bytes\n", g_transfer_index_bytes);
if(g_transfer_data_bytes > 1024*1024)
Log(" Data : %11u : %7.2fM bytes\n", g_transfer_data_bytes, g_transfer_data_bytes/(1024.0f * 1024.0f));
else if(g_transfer_data_bytes > 1024)
Log(" Data : %11u : %7.2fk bytes\n", g_transfer_data_bytes, g_transfer_data_bytes/1024.0f);
else
Log(" Indices : %11u bytes\n", g_transfer_data_bytes);
}
#endif